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We demonstrate that graphs embedded on surfaces are a powerful 

and practical tool to generate, characterize and simulate networks with 

[~^. ■ a broad range of properties. Remarkably, the study of topologically em- 

{^J ' bedded graphs is non-restrictive because any network can be embedded 

on a surface with sufficiently high genus. The local properties of the 
network are affected by the surface genus which, for example, produces 
significant changes in the degree distribution and in the clustering coeffi- 
cient. The global properties of the graph are also strongly affected by the 
f\ ' surface genus which is constraining the degree of interwoveness, changing 

*~i , the scaling properties from large-world-kind (small genus) to small- and 

ultra-small- world-kind (large genus) . Two elementary moves allow the ex- 
ploration of all networks embeddable on a given surface and naturally in- 
troduce a tool to develop a statistical mechanics description. Within such 
a framework, we study the properties of topologically-embedded graphs 
at high and low 'temperatures' observing the formation of increasingly 
regular structures by cooling the system. We show that the cooling dy- 
namics is strongly affected by the surface genus with the manifestation of 
a glassy-like freezing transitions occurring when the amount of topological 
disorder is low. 
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1 Introduction 

In natural and artificial systems there exists a very broad variety of networks; 
indeed, networking is a crucial feature in information technologies, it is a vital 
skill in social behaviour and -more generally- it is at the base of the emergence 
of some of the fundamental properties of complex systems. Networks possess 
a wide range of properties and their structure can assume different forms de- 
pending on their function, their construction rules and their evolution dynamics. 
From a general perspective, the structural properties of a network can be divided 
into two main categories: (i) the local structure, concerning small portions of 
the graph which may vary from place to place and, typically, they are analyzed 
statistically (a well-known example is the degree distribution [T] : distribution of 
the number of edges per vertex); (ii) the global structure, which concerns prop- 
erties that involve the entire organization of the graph (a well-known example 
is the diameter [2d]: the longest shortest path between any two vertices). It is 
understood that local properties and global properties are related, for instance, 
it has been shown that the degree distribution is affecting the diameter and, 
in particular, the presence of a few highly connected hub- vertices can reduce 
significantly the overall diameter [HE]- However, the relation local/global is in 
general mediated through the hierarchical organization of the network and can 
result in non-trivial relations. By their nature, local properties are easier to be 
measured and therefore, so far, they have attracted most of the attention in the 
literature. 

In this paper we show that by considering networks embedded on surfaces 
we can control an important global measure of complexity which is associated 
with the number of cycles: the network interwovenness. It is quite intuitive 
that an increased interwovenness must be related with an increase of the com- 
plexity of the network structure, however, due to its global nature, a general 
quantification of this quantity is a very challenging task. On the other hand, 
by considering the embedding of a network on a surface the interwovenness can 
be directly associated with the surface genus which is a non-negative integer 
number counting the number of handles in the surface [3]. Let us recall that 
a network topologically embedded on a surface has the property that it can be 
drawn on the surface without edge-crossings (edges do not have to be straight). 
The absence of edge-crossing is obviously a limitation on the degree of interwo- 
venness of the network and it is therefore a constraint on its overall complexity. 
A sphere has no handles and genus g = 0, a torus has one handle and g = 1, a 
double torus has g = 2, etc.. Intuitively, we can look at a handle in a surface 
as a 'short-cut' that connects two distant parts. For instance, by joining the 
north and south poles of a sphere and pinching them together one can trans- 
form the sphere into a tours passing from g = to g = 1. It is clear that these 



short-cuts can be also used by the embedded network that in this way can link 
otherwise distant vertices. In this paper we consider only simple graphs, where 
no more than one edge can directly connect two vertices and there are no loops 
connecting a vertex to itself. If we take an orientable surface and we place n 
vertices on it, we can then connect with edges couples of vertices up to a point 
when no further edges can be inserted without generating edge-crossings. We 
call this graph maximal embedded graph [4 and, with the exemption of some 
special cases, it is a triangulation of the surface containing 3n + 6(g — 1) edges. 
By increasing the genus we can insert an increasingly larger number of edges in 
the maximal embedded graph. The genus is therefore a measure of the graph 
interwovenness. Any other graph with n vertices embedded on the same surface 
must be a subgraph of a maximal embedded graph (indeed, either it is a max- 
imal embedded graph, or edges can be added to make it a maximal embedded 
graph). Furthermore, it has been proved that any graph can be embedded on 
a orientable surface with sufficiently large genus [5], therefore, by considering 
maximal embedded graphs we are embracing the whole family of all possible 
graphs. 

One of the advantages of considering maximal embedded graphs is that there 
is a simple constructive way to build and modify them. In general, there exist a 
large number of models and construction methods that allow to build networks 
with desired properties and structures [TJ[2] • Ideally, it would be desirable to be 
able to consistently generate networks with controlled and tunable properties 
both at local and global levels. Such networks must be able to evolve and adapt 
following simple mechanisms and, eventually, allowing for a statistical mechanics 
kind of approach to be implemented. Maximal embedded graphs can be easily 
generated by starting from a seed embedded structure and then they can evolve 
by means of two elementary moves, known as Tl and T2 6-8 , which involve 
only local changes and do not modify the embedding (see FigfT]). 

The paper is organized as follows: in section [2] we introduce the idea of em- 
bedding networks on surfaces, we discuss the properties that embedded networks 
must obey, we describe the elementary moves Tl and T2, we define a simple 
energy function and we put the basis for a statistical mechanics description of 
these systems. In section [3] we perform an extensive numerical study of the 
properties of maximal embedded networks discussing the effect of the surface 
genus and of the temperature on the structural and dynamical characteristics 
of these complex networks. In section|4l conclusions and perspectives are given. 

2 Maximal Embedded Graphs 

2.1 Topological properties 

As mentioned earlier, a maximal embedded graph on a surface of genus g is a 
triangulation of the surface. If n is the number of vertices, u the number of 
edges and t the number of triangles, then the Euler's polyhedron formula [9 




Figure 1: Exemplification of the two elementary moves which allow to explore 
all possible maximal embedded graphs on a given surface, (a) the Tf move 
which consists in the switching of one edge between four vertices (a, b, c, and 
d); (b) the T2 move, consisting in the insertion of a vertex (d) inside an existing 
triangle. 



provides us with a very simpe relation between these numbers 

n-u + t = 2(1 - g) . 



(1) 



Furthermore, we have that, in average, each vertex has (fc) incident edges and 
each edge has two vertices at its extremities (recall that we consider simple 
graphs that do not allow multiple edges between couples of vertices) , yielding 
to 

(k)n = 2u . (2) 

Moreover, each triangle has three edges and each edge is in between two trian- 
gles, however some triangles may be self-neighbours meaning that the edge may 
have the same triangle in both sides, implying 



3t = 6u 



(3) 



where 9 € (0, 2] with 9 = 2 associated to the case when triangles cannot be 
self neighbours and the limit 9 — » associated to the case when all triangles 
are self-neighbours. By combining the previous two expressions with the Euler 
formula, EqUJ we get 



2 + 4^-^<(fc)<6 + 12^ — - 



(4) 



where the upper bound is achieved only when triangles cannot be self-neighbours. 
Incidentally, EqU provides us with a lower bound for the genus required 
to embed any network. Indeed, by construction, any network embedded on 
an orientable surface of genus g must be a subgraph of a maximal embedded 
graph and therefore it must have an average degree smaller or equal than it. 



Therefore, if k is the average degree of an arbitrary network (not necessarily a 
triangulation) , then from Eq|4]we have 

fc-6 
9 > 1 + -jg-n • (5) 

Although, this bound is achieved by maximally embedded graphs with no self- 
neighbouring triangles, it is in general a rather loose bound for complex net- 
works. We can, for instance, note that for a large sparse network with in average 
less than six edges per vertex, the right-hand side of Eqj5]becomes negative, but 
g must be a non-negative integer and therefore the bound becomes ineffective. 
It must be noted that this bound can be also applied locally to any sub-graph. 
This may give more restrictive bounds especially when highly compact sub- 
structures are present. For instance, if we have a q-clique (a configuration of q 
vertices each-one connected with all the others), then for this sub-structure we 
can substitute k = q — 1 in EqJS] obtaining g > {q — 2>){q — 4)/12, which is the 
same bound proved by Ringel and Young in [5j[l0] for the solution of the map- 
coloring problem. For instance, we can see that a tetrahedral clique {q — 4) can 
be embedded on a sphere (g > 0), but a 5-clique (K 5 , the complete graph with 
five vertices) requires already a larger genus g > 1. This is in agreement with 
the Kuratowski theorem [llj which proves that planar graphs cannot contain 
K$ or K33 (the complete bipartite graph with six vertices connected three by 
three) as minor. 

2.2 Elementary moves 

Maximal embedded graphs can be built by starting from a seed-graph and then 
by letting it evolve through simple elementary moves, called Tl and T2 [MS]- 
The first move consists in the switching of an edge in a local configuration of 
four vertices in which two second neighbour vertices become directly connected 
and vice-versa two first neighbour vertices become second neighbours. This is 
shown in FigQTa). New vertices can be added to the graph by means of the 
second move, T2, which consists in the insertion of a vertex within an existing 
triangle generating in this way three new triangles, as shown in FigJIJb) [6]. 
Similarly, vertices can be removed from the system by applying an inverse T2 
move eliminating three triangles included inside a 3-clique. These are local 
moves that do not change the global embedding of the graph. For a given 
surface, by means of Tl and T2 moves, it is easy to build and explore the entire 
class of maximal embedded graphs [5][7] . Interestingly, graph generated by T2 
moves only are known in the network-litereature under the name of Apollonian 
networks [12J . 

2.3 Energy 

We have now the tools to develop a statistical physics framework for maximal 
embedded graphs. To this purpose it is convenient to introduce a topological 



energy which can guide us in the exploration of configurations with different 
degrees of disorder. A simple form for the energy is: 

n 

E = Y / (h-(k)) 2 , (6) 

i=i 

which associates the ground state (E = 0) to the configuration with lowest 
topological disorder where all vertices have the same degree k = (k) (if (k) 
is integer). This energy associates a 'cost' to the deviations from the average 
degree (fc) and it is therefore a distance measure (the square of an Euclidean 
distance) from an ideal regular network with degree (k). 

Our approach is to use a combined action of Tl and T2 moves to build a 
maximal embedded graph with a given number of vertices n on a surface with 
given genus g . Once the graph is built we explore the achievable configurations 
by means of Tl moves only. 

After a Tl move, two vertices (i.e. k a and k c in FigOJa)) acquire a new edge 
and the other two vertices (i.e. kb and kd in FigOJa)) loose one edge leaving the 
overall (k) unchanged. As a consequence of such a move the energy change is: 

AE = (k a + 1 - <fc» 2 + (*c + 1 - (k)) 2 + (h-l- (k)f + {k d - 1 - (k)f - 

- (k a - (k)) 2 - (k c - (k)f - (k b - (k)) 2 - (k d - {k)f , (7) 

which is 

AE = 2(k a + k c -k b -k d )+4 . (8) 

Intriguingly, the energy variation as consequence of a Tl move does not depend 
on (k) (i.e. it is independent on the embedding) and it is linearly dependent on 
the degrees of the four vertices involved in the Tl move. One can immediately 
see that the Tl move decreases the system-energy iff k a + k c < kb + k d ; mean- 
ing that the energy is decreased if the nodes with smaller degree acquire new 
connections at the expenses of the nodes with larger degree. We are therefore 
describing a local averaging process where inequalities in the degrees are grad- 
ually redistributed lowering in this way the total energy and driving the system 
towards the ground state. 

2.4 Dynamics 

A 'temperature' /3 _1 can be introduced and a statistical mechanics description 
can be implemented. To this end, we adopt a Glauber-Kawasaki type of dy- 
namics where a Tl move is performed accordingly with the probability |13j 

n(fc„, h, kc, kd)= 1 + ex p (/3Ajg) (i - **»,3)(1 - fod.3)(l - Kc) , (9) 

where the Kronecker deltas enforce a 'proper' triangulation. Specifically, the 
first two are preventing the vertex degree to become smaller than 3 and the 
third avoids the formations of loops where edges connect a vertex to itself. 
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Figure 2: Average energy per vertex (E) jn for the cooled states after 10 4 x n 
attempted Tl moves at a given temperature /3 _1 for maximal embedded graphs 
with n = 20000. Symbols: □ g = n + 1; Q) g = 0. The dashed lines are the 
mean-field predictions for the two genera. 



2.5 Mean- field solution 

The maximal embedded graphs that we are studying here have a very simple 
energy that depends only on the degree of the nodes. Nevertheless, interac- 
tions between vertices arise unavoidably from the network topology and its 
constraints, making any analytical study extremely challenging [14H19J . How- 
ever, at low enough temperatures, the energy should dominate over entropy and 
it may be possible to describe the equilibrium state of the system by assuming 
all vertices as independent and reducing the topological correlations to a sim- 
ple constraint on the average degree. In such a mean field approximation the 
partition function factorizes and it can be written as: 




exp 



(-/3(fc-(fc)) 2 



+ Xk) 



(10) 



where the coefficient A is a Lagrange multiplier associated to the average degree 
and it is (implicitly) given by the condition: 



d 
(k) = ttt In Z n . 

The average energy at equilibrium is given instead by 

{E) = -^lnZ n , 



(11) 



(12) 



which can be calculated numerically. An analytical expression can be obtained 
in the limit of low temperatures revealing that the energy should decrease lin- 
early with the temperature: (E) jn ~ /3 _1 /2. In the opposite limit of high 
temperatures, we have instead that the energy is expected to tend to a con- 
stant: (E) /n~ ((k) -3) 2 . 

3 Numerical results 

We have numerically generated maximal embedded graphs containing up to 
n = 20000 vertices embedded on surfaces with different genera ranging between 
g = to g = 2n + 1. 

3.1 The random states 

Let us start from high temperatures where the transition probability, Eqj9l 
becomes independent on the degrees of the four nodes (except for the forbidden 
moves) . After a large number of moves the resulting equilibrium configurations 
become statistically stable and we shall refer to these configurations as the 
random states. For g = some properties of this state are analytically known 
[T4HT9] . For instance, the degree distribution is: p(k) = 16(3/16) fe (fc — 2) (2k — 
2)!/[fc!(fc — 1)!], which, in the tail region, is well described by an exponential 
law. This is in a satisfactory agreement with the mean field description that 
implies p(k) — ((k) — 3) fe ~ 3 ((fc) — 2) 2 ~ fc . However, when the genus increases, 
the distribution deviates from the exponential acquiring a power law behaviour 
in the region k < (k) and simultaneously assuming a faster than exponential 
decay in the region of large degrees k > (k). We previously reported in [20] that 
the degree distribution for arbitrary genus appears to be well described by the 
following functional form p(k) oc fc~"cxp(— Afc 7 ) where the parameters depend 
on the genus with a passing from zero at g — to about 1 at <? = n + 1, A 
changing from about 0.3 at g = to about 0.03 at g = n + 1 and 7 passing from 
about 1 at g = to about 2 at g = n + 1 4 . These changes of the parameters 
reflect the change in the degree distribution from exponential to power-law with 
a faster than exponential cut-off at large k. A detailed analysis of the properties 
of this state is reported in [4] . 

3.2 Equilibrium at finite temperatures 

In FigJ5] we report the average energy measured after 10 4 x n attempted Tl 
moves at a given temperature /3" 1 on maximal embedded graphs with n = 20000 
vertices. These values are averages over 10 cooling loops and the data refer to 
two different embeddings respectively with 5 = and g — n + 1. The dashed 
lines are the mean-field predictions from Eq[12] (calculated in the limit moo 
and by substituting the sum with an integral). One can note that there is a 
range of temperatures in which the average energy follows very well the mean 
field prediction. This is the regime where the network is at thermal equilibrium 



and the quadratic energy is dominating the degree distribution which is there- 
fore well approximated by a normal distribution. We can observe that at high 
temperatures (/3 S> 1), the average energy saturates to a plateau following a 
qualitatively similar behaviour to the mean field prediction but with significant 
quantitative deviations. Indeed, at high temperatures, the statistical properties 
of the network are dominated by the entropic term which is associated to the 
number of triangulations that can be built on a surface of a given genus [14] . 
On the other hand, the mean field equations take into account the topological 
properties of the network only by constraining the average degree, and this turns 
out to be not sufficient for quantitative estimates. 

3.3 The freezing dynamics 

Starting from high temperatures (/3 _1 3> 1) we can observe from FigJ5]that the 
residual energy decreases with the temperature indicating that -as expected- 
the states at lower temperatures are less disordered. As mentioned earlier, we 
observe that at intermadiate temperatures (/3 _1 ~ 1) the system follows the 
mean field prediction, however, around j3~ l ~ 0.5 a change in the behaviour 
occurs. At low temperatures, the residual energy does not longer decrease with 
the temperature approaching instead a constant plateau. This is the signature 
of a 'freezing' phenomena [13l[211|24] typical of glasses [T3ll25ll26] . 

3.4 The frozen states 

At low temperatures (/3 _1 <C 1) the maximal embedded graph should approach 
the ground state where all vertices have degree equal to (k) (if it is an integer 
number). However, we observe that the system cannot reach the ground state in 
any finite time. This freezing phenomena was first observed in studies performed 
on planar triangulations with periodic boundary conditions (g = 1) [13,21-24 
and here it is retrieved in the general case of maximal embedded graphs on 
surfaces of arbitrary genus. Interestingly, we observe that the complexity of 
the embedding surface (e.g. its genus g) is affecting the freezing dynamics 
and consequently the properties of the asymptotic states @|. In particular, 
triangulations with large genus start from random states with higher energies 
with respect to the low genus counterparts, but they cool faster and, within the 
same number of attempted Tl moves, they can reach frozen states with lower 
residual energy. When the average degree is an integer, the degree distribution of 
the frozen states is characterized by a large fraction of vertices at a degree equal 
to the average and by two small fractions of 'defective' vertices respectively with 
degree one above and one below the average (see [4]). A similar distribution 
is observed when the average degree is not an integer but it is still close to 
a natural number. On the other hand, a different statistics and a different 
dynamics are observed when (k) is near to an half-integer, an effect due to 
topological frustration. 
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Figure 3: between the energy at different times (E) and the expected energy of 
the ground state {Eoo) divided by the average energy of the random state (Eq). 
Symbols: O g = 1, (k) = 6; U\ g = 10000, (k) ~ 12; o g = 20001, (k) = 18; > 
g = 834, (k) ~ 6.5; V g = 12501, (k) = 7.5; A g = 17501, (A;) = 16.5; 

3.5 The topological frustration 

When (k) is in the middle between the two integers, k\ and k%, then the energy 
is minimized by a state with half of the vertices with degree equal to k\ and 
the other half with degree equal to ki . This naturally introduces frustration in 
the system by preventing the formation of a regular ground state. In particular, 
the ground state energy passes from E^ — for (fc) integer to E^ — 0.25 n 
for (k) half integer. Intriguingly, we observe that this also strongly affects the 
dynamics. In Fig[3] we compare the cooling dynamics in three networks with 
nearly integer average degrees (namely g = 1, 10000, 20001 which correspond 
respectively to (k) ~ 6, 12, 18) and in three other networks with nearly half- 
integer average degrees (namely g = 834, g — 2501, g = 17501, corresponding to 
(k) ~ 6.5, 7.5, 16.5). In all cases, the system is prepared in the random state at 
infinite temperature (/3 = 0) and then the temperature is quenched to zero. The 
figure reports the difference between the values of the energy during the cooling 
dynamics after quenching and the expected energy of the ground state (Eoo) 
divided by the average energy of the random state (E ). In this way, all plots 
start from 1 and they must approach if the system is reaching equilibrium. 
We observe that the unfrustrated systems with integer (k) freeze to a plateau 
energy which decreases with the genus but it is well above Eq = for all the 
studied cases. We instead observe that the frustrated systems with half-integer 
(k) exhibit a fast dynamics to E^ — 0.25 n which is reached in finite times with 
exponentially fast relaxation times. 
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Figure 4: Variation of the graph-diameter vs. the genus. Symbols: A random 
states; V frozen states. 



3.6 Shortcuts 

Let us now investigate the effect of the surface-genus on the global properties of 
maximal embedded graphs. As mentioned earlier, the presence of an handle in 
the surface provides a possible pathway for 'shortcuts' between two distant parts 
of the network that otherwise could not be directly connected. It is intuitive to 
understand that by increasing the genus one also increase the number of short- 
cuts that can be inserted. Consequently the diameter of the graph (maximum 
distance between any two nodes) must proportionally decrease when the genus 
increase. However, we must stress that we are investigating a stochastic system 
that is spontaneously evolving through random weighted moves and therefore 
the relation between network-diameter and surface-genus is not straightforward. 
In Fig [J] we report the graph diameters for samples with n = 20000 vertices em- 
bedded on surfaces with different genus. As expected, we can see that the 
diameter decreases with the genus. We observe that this happens both for net- 
works in the random and in the frozen states. At low genus the random state has 
a smaller diameter than the frozen state, a feature that may be associated with 
the formation of hubs in the random configuration. Intriguingly, we can see that 
at high genus the relation is reversed and frozen states have smaller diameters 
than random states, this may be associated to the formation of branches in the 
random state. 
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Figure 5: Variation of the graph-diameter vs. the number of vertices. Symbols: 
A random states; V frozen states. 
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Figure 6: Variation of the graph-diameter vs. the number of vertices. Symbols: 
A random states; V frozen states. 
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3.7 Scaling of the graph diameter with network size 

We have also investigated the effect of the surface genus on the scaling properties 
of maximal embedded graphs. Networks can have different scaling properties 
which can be quantified by estimating the dimension [27] . For instance, a 
triangular lattice on a two-dimesional Euclidean surface has a dimension equal to 
two. This can be measured by investigating the scaling of the network-diameter 
with the number of vertices in the network which, in this case, must scale with 
the square of the diameter and therefore we have the relation: diameter oc n 1 ' 2 . 
More generally, if the relation between the diameter and the size is a power law, 
then one can define a dimension d from: diameter oc n}' d . 

In Figf5]the graph-diameter is plotted vs. the number of vertices for maximal 
embedded graphs with low genera {g = 0) and high genera (g = n + 1) both for 
random and frozen states. As clearly visible from the plots, the genus strongly 
affects the value and the scaling of the diameters. 

Incidentally, the dimension of maximal embedded graphs on spherical sur- 
faces (g — 0) has been studied in the context of two-dimesional quantum gravity 
where the (Euclidean) space-time is associated to a dynamical triangulation with 
all edge-lenghts constant but with variable vertex degree which is associated to 
the local space-time curvature [2"5H3"T] . Studies in this context have revealed that 
the dimension must be equal to 4 [32 j ; a result comforted by an analytical calcu- 
lation via a transfer matrix formalism 33 . In the present study we observe that 
the diameter follows rather well the functional form: diameter ~ cin 1 ' 4 + c%. 
However, by looking carefully at Figj5] one may note that the power law be- 
haviour is satisfactory for g = but, instead, for high genus cases, the growth is 
slower than any power law. This is made evident in Fig|6] where the diameter 
is plotted as function of ln(n) and ln(ln(n)). We can see clearly in FigfBJa) that 
the diameter in low genus networks grows faster than ln(n) indicating a power 
law scaling. On the other hand, we can see from FigOJb) that the diameter in 
high genus networks does not grow any faster than ln(n) and the growth may 
even be consistent with a linear increase with ln(ln(n)) suggesting that these 
networks could be ultra-small worlds. Ultra-small world properties have been 
observed in scale-free networks |34j but here -instead- this property is revealed 
also in the frozen state which are almost regular graphs with very narrow degree 
distributions. It must be stressed that it is beyond the purpose of the present 
work to establish the exact scaling laws, indeed, for this purpose, a larger range 
of sizes must be explored. However, what we have observed gives enough qual- 
itative evidences to establish that the surface genus strongly affects the global 
properties of maximal embedded graphs changing the scaling of the diameter 
with the number of vertices passing from a large-world-kind of behaviour to a 
small-, and possible even ultrasmall-world behaviour. 

3.8 Clustering coefficient 

Another important and widely used graph-measure is the clustering coefficient 
which quantifies the level of local interrelations between a set of neighbouring 
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Figure 7: Variation of the average clustering coefficient vs. the genus. Symbols: 
A random states; V frozen states. The line is the prediction from c = 2(5 + 
12(g- f)/?i) _1 (see text). 
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vertices. We measured the clustering coefficient of each vertex by taking the set 
of its neighbouring vertices and counting the number of edges between them and 
then dividing this number by the maximum possible number of edges between 
them (which is k(k — l)/2 for a vertex with degree k). In FigJT] we plot the 
average clustering coefficient for both random and frozen states. We can observe 
that maximal embedded graphs have finite clustering coefficient which decreases 
with the increasing genus. This is related to the fact that the average number of 
neighbours increases with the genus (EqHJ) . Indeed, in a triangulation, a vertex 
with degree k has at least k edges between its neighbours. This means that, for 
a vertex with degree k the clustering coefficient must have the following lower 
bond: 

*> * i(£lj = F^T • < 13 » 

which is a decreasing function of k. The equality in Eqll3l is achieved when 
the local structure surrounding the vertex with degree k is a polyhedron with 
k vertices. If the degree distribution is narrowly distributed around (k), then 

we can substitute ( ^-j \ ~ (k)-i an< ^ therefore from Eqs. [4] and [13] we can 
derive a lower bound on the average clustering coefficient: c = / k \_ 1 = 2(5 + 
12 (g — l)/n) _1 . This lower bond is plotted with a continuous line in FigfTl As 
one can see, the line overlaps well with the results for the frozen states. This 
is intuitively correct because the frozen states are almost regular graphs with 
planar local configurations. 



4 Conclusions 

We have shown that by considering networks embedded on surfaces one can 
develop a powerful formalism that allows to study any kind of network within 
a unified statistical mechanics approach. The surface genus is a very impor- 
tant characteristic that strongly influences both local and global properties of 
the embedded graph. There are two simple, local elementary moves, Tl and 
T2, that allow to explore the entire set of maximally embedded graph on a 
given surface. By means of these moves one can develop a statistical mechanics 
approach which can be used to investigate networks with different degrees of 
disorder but with constrained complexity. At high temperatures, the network 
is disordered. At low genera the equilibrium random state reveals an exponen- 
tial degree distribution. Instead, at high genera, the equilibrium distribution 
becomes a power law distribution with faster than exponential cut-off. Con- 
versely, at low temperatures, the network has an almost regular degree with 
only few 'defective' vertices. The cooling dynamics incurs in a slowing down 
when the amount of disorder in the structure becomes low and the ground state 
cannot be reached in a finite time. Intriguingly, we observe that a mechanism 
of geometrical frustration, associated with fractional average degree, can rein- 
troduce fast cooling dynamics showing that the occurrence of this glassy phase 
is related to the amount of disorder. We show that a mean field model can 
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describe correctly some of the network properties especially in the range of fi- 
nite temperatures. From a global structural perspective, the surface genus is 
constraining the degree of interwovenness of the network and it is affecting the 
graph-diameter. We found that networks developing on surfaces with small gen- 
era have large- world properties with the diameter scaling as a power law of the 
number of vertices, instead high genera surfaces lead to much smaller networks 
that might be small-world or even ultra-small- worlds. Interestingly, this scaling 
is observed in networks prepared at high temperatures and possessing broad 
degree distributions as well as in networks cooled at low temperature with very 
regular structures. 

Future studies will focus on the application of these networks to information 
filtering by constructing maximal embedded graphs from a similarity measure 
in analogy with what already done with the PMFG construction [35] that in 
the present context corresponds to the g — case. 
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